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OSCILLATIONS  OP  A  PHYSICAL  PENDULUM  HAVING  CAVITIES 
FILLED  WITH  A  VISCOUS  FLUID 


P.  S.  Krasnoshchekoy 
(Moscow) 

An  asymptotic  method  has  been  proposed  [1]  for 
investigating  the  nonstationary  motions  of  viscous 
incompressible  fluids  with  large  Reynolds  numbers 
(NRg),  which  arise  during  the  oscillation  of 

various  solids  which  either  contain  a  fluid  or 
are  emersed  in  one,  and  also  during  the  oscilla¬ 
tion  of  liquid  volumes  having  a  free  surface. 

In  this  article  the  idea  of  this  method  is  applied 
to  the  investigation  of  small  oscillations  of  a 
physical  pendulum  having  cavities  completely  filled 
with  an  incompressible  viscous  fluid. 


1.  The  pressure  of  a  fluid  contained  in  the  cavity  of  an  oscil¬ 
lating  pendulum  (Fig.  1)  Is  described  by  the  Navler-Stokes  equations 
and  the  continuity  equation: 

W  +  <v'  •  v'>  v'  =  v  V  -  v  (V'  x  O') 

divV'  =  0  '  (1-1) 

(a'  =  v'x  v',  <p'  =  -£-*/) 

Here  U  is  the  potential  of  the  mass  forces  acting  on  the  fluid. 

At  the  edge  of  the  cavity  the  conditions  of  adhesion  of  fluid 
particles  to  the  walls  of  the  cavity  should  be  satisfied.  This  gives 
the  following  boundary  conditions: 
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=  —  j/'Q,  v'  =  z'0, 


=  0 


(1.2) 


Here  8  is  the  deflection  angle  of  the  pendulum  from  the  equilib¬ 
rium  position  while  u1,  vr,  w’,  are  the  components  of  the  velocity 
vector  y . 

We  assume  that 


y  =  v"  +  Yq 

where  Yq  is  the  velocity  vector  of  the  center  of  mass  of  the  cavity, 
and  we  convert  to  a  new  system  of  coordinates  (x",  y",  z" )  with  its 
origin  at  the  center  of  mass  of  the  cavity  and  with  axes  parallel 
to  the  axes  of  the  stationary  system  of  coordinates  (x1,  y1,  z1)* 

In  the  new  system  of  coordinates  we  will  have 

3F+ (V'.V')V' =  VV-v(V'xQ').  '  divV*-0.,  (1>5) 

Here 

Q'»V'xV', 

At  the  boundary  of  the  cavity: 

.  -  y- 6‘,  v'  =  z'0,  w"  =  0  (1.4) 

We  will  refer  all  quantities  to  a  characteristic  scale.  Let  us 


assume: 


t'  =  37,  .  z'  -  Rx,  y’  =  Ry , 


=  V'  =  3-Ry 

0  =  ad,  JVRe  =  £ 


(1.5) 


where  T  is  the  characteristic  period  of  one  oscillation,  R  is  the 
characteristic  dimension  of  the  cavity,  a  is  the  characteristic  ampli¬ 
tude,  and  NRe  is  the  Reynolds  number. 

Equations  (1.5)  and  Boundary  Conditions  (1.4)  are  written  in 
independent  variables  in  the  form: 


^  +  a(V.V)V  =  vcp  — ^(VxQ). 

R® 


div  V 


(1.6) 
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At  the  boundary  of  the  cavity: 

“  =  —  y&,  V  =  2T&,  W  .=  o  (1*7) 

In  the  system  of  equations  (1.6)  the  characteristic  amplitude  a 
is  a  factor  for  nonlinear  terms. 

We  will  consider  oscillations  with  a  small  amplitude  and  later 
linearize  Eqs.  (1.6),  disregarding  magnitudes  of  the  order  a. 

The  linearized  system  of  equations  appears  as  follows: 

~  =  Vip— J-(Vxft),  div  V  =  o  (1*®) 

'Re 

The  object  of  this  article  is  to  investigate  those  forms  of 
oscillations  for  which  the  solution  may  be  represented  in  the  form: 

*''■)  =  «•'.  V  =  (X,  y,  Z)  (1.9) 

Setting  q>  =  ee-'O  (x,  y,  z),  =  ce”'^  (x,  y,  z),  we  have  for  such  motions: 

=  VO-^fVxip),  div  U  =  0  (1.10) 

At  the  boundary  of  the  cavity: 

Ux  =  -  yo,  Uv  =  Xc,  Ur  =  0  (1.11) 

Relationship  (1.10)  indicates  that  the  vector  TJ  is  the  stun  of  the 
potential  and  solenoidal  vectors.  Such  a  presentation  makes  it  possible 
to  separate  the  equations,  having  obtained  a  separate  equation  for 
each  unknown  function. 

Actually,  having  taken  the  div  of  both  sides  of  the  first  equation 
of  System  (1.10)  we  obtain: 

AO  =  o  (1.12) 

i.e.,  a  function  of  <1>  which  is  harmonic  in  a  volume  filled  with  a  fluid. 
Having  taken  the  cure  of  both  sides  of  the  same  equation  we  .have: 

(i.w) 

Thus  each  of  the  unknown  functions  $>,  Y  satisfies  a  separate 
equation. 
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At  the  boundary  of  the  cavity  these  functions  are  linked  by  the 


following  boundary  conditions: 


a® 

1  1 

r  _ 

ox  - 

l  Sy 

dz  j 

3® 

1  i 
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=  xo1 
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Here  Y  .  Y  ,  Y  are  components  of  vector  Y, 

X  jT  £ 

A  fourth  boundary  condition  determined  at  the  limit  of  the 
projection  of  vector  'S  onto  the  normal  to  the  surface  of  the  cavity 
will  be  introduced  later  in  a  specially  selected  curvilinear  system 
of  coordinates. 

Let  us  consider  cavities  which  appear  as  figures  of  revolution 
whose  sixes  are  perpendicular  to  the  plane  of  oscillations  of  the  body. 

We  will  introduce  a  curvilinear  system  of  coordinates  associated 
with  the  surface  of  the  cavity.  In  the  case  under  consideration  it 
is  convenient  to  introduce  the  following  coordinates:  n  —  distance 
along  the  inner  normal  to  the  surface  of  the  cavity  taken  from  the 
surface  Inward,  <3  —  angle  which  determines  the  position  of  the  meridian 
plane,  and  P  —  length  of  arc  along  the  meridian  (Fig.  2).  The  varia¬ 
bles  in  the  (x,y,z)  coordinate  system  are  associated  with  the  variables 
n,  a,  P  in  the  curvilinear  coordinate  system  by  the  following  relation¬ 
ships: 

a 

*  =  -  [r0  +  $  sin  r  (P)  -  n  cos  r  (P)J  sin  a  =  -  l  r  (P)  -  n  cos  T  (P)l  <* 

0 

0 

y  =  [r0  -f  $  sin-T  (P)  -  n  cos  r  (P)J  cos  a  -  [ r  (P)  -  n  cos  r  (P)l  cos  0 

0 

0 

z  =  ^  cos  I  (M  $  +  n  sin  T  (P) 
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Here,  rQ  =  radius  of  a  circle  obtained  in  the  cross-section  of  the 
v  cavity  by  a  plane  perpendicular  to  the  axis  of  symmetry 
and  passing  through  the  center  of  mass  of  the  cavity; 

r(0)  =  radius  of  the  circle  formed  by  points  on  the  surface  of 
the  cavity  having  the  coordinate  0; 

7(0)  =  angle  between  the  tangent  to  the  meridian  at  the  point  with 
coordinate  0  and  the  positive  direction  of  the  z-axis. 

We  will  write  Eq.  (1.13)  and  the  boundary  conditions  in  the  new  coor¬ 
dinate  system.  For  simplicity  we  will  use  the  notation: 

v*  =  rv  -£ 


Then  we  have  in  the  curvilinear  coordinate  system: 


+ 


cw,  _ _ 1 _ ro-/  l±*£M  dVl) 

■  ^Re  lr  (3)  —  n  COS  T  (3))  [1  +  nt'  (3)]Laa  \r  (3)  —  n  cos  y  (8)  da  J  ' 

+  m  ^) + r. +<«-■>»’  r  MUi  +  «m  5)] 

(f  (»-%“' 


(i  =  1,  2,  3) 


At  the  edge  of  the  cavity  for  n  =  0: 
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! 

dn 
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(1.14) 


(1.15) 

( 1 . 16 ) 

(1.17) 


Let  us  derive  the  fourth  boundary  condition.  The  normal  compo¬ 
nent  ¥n  of  the  vortex  vector  Y  is  defined  in  terms  of  the  tangential 
components  UQ  and  Up  of  the  velocity  vector  U  as  follows: 

' iff _ 1 _ fd  U  +  "f  (3))  ua  d  (r  (3)  -  n  cos  T  (3))  CUa 

"  (<•  (P)  —  «  cos  r  (.3)1  (1  +  ny'  (3)]  [  5a  :  53 - — J 


At  the  edge  of  the  cavity  the  vector  U  is  known;  consequently, 
its  components  UQ  and  Up  are  known  for  n  =  0.  Differentiation  with 
respect  to  a  and  0  for  n  =  const  is  possible,  and  therefore  calcula¬ 
tion  of  the  magnitude  of  Yn  for  n  =  0  presents  no  particular  difficulty 
After  simple  calculations  we  obtain: 
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2  sin  r  (p)  o 


=  - 


(1.18) 


This  relationship  closes  the  system  of  boundary  conditions  for 
the  unknown  functions  Y^. 

The  quantities  YQ,  Y^,  Y n  in  Relationships  (1.15)-(i.l8),  are 
components  of  vector  Y  in  the  curvilinear  coordinate  system.  They 
are  connected  with  the  componenets  Y^  by  the  relationships: 

=  —  Y,  cos  a  —  sin  a 

V>  =  -  sin  T  (P)  sin  a  +  Y,  sin  T  (p)  cos  a  +  Y,  cos  r  (p)  (l.  19) 

Yn  =  Y,  cos  r  (P)  sin  a  -  cos  r  (P)  cos  a  +  sin  r  (P) 


We  will  assume  that  the  parameters  of  the  pendulum  guarantee  a 

(1.20) 


sufficiently  large  Reynolds  number.  We  will  assume 

i 


7T~  =  r2 

nRe 


where  e  is  a  small  dimensionless  parameter. 

The  idea  of  constructing  a  solution  for  large  Reynolds  numbers 
which  was  presented  by  Moiseyev  [l]  is  analogous  to  the  idea  of 
constructing  a  boundary  layer.  It  is  assumed  that  the  vortices  which 
exist  in  the  oscillating  fluid  contained  within  the  pendulum  cavity 
are  concentrated  mainly  in  a  thin  layer  at  the  walls  of  the  cavity. 
This  in  turn  makes  it  possible  to  consider  the  derivatives  of  the 
components  of  the  vector  Y  along  the  normal  to  the  surface  of  the 
cavity  to  be  significantly  greater  than  in  the  tangential  directions. 
Let  us  "expand"  the  independent  variable  n.  We  set 

n  =  et) 

Solution  of  the  problem  formulated  will  be  sought  in  the  form  of 
series  to  powers  of  the  small  parameter  e: 

*  -  ®,  +  e®,  +  .  .  .,  Y,=±y#t+  Ylt  +  ...  (i-1.2.3)  (1.21) 

Substituting  Series  (1.21)  Into  the  equations  and  boundary 
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conditions  and  equating  to  zero  the  sum  of  the  coefficients  for  the 
zero  power  of  the  parameter  e,  we  obtain  the  following  problem  for 
the  determination  of  the  functions  3>q  and  : 

A<X>0  =  0,  (7^  =  ^  <i  =  1, 2. 3)  (1.22) 


At  the  boundary  of  the  cavity,  for  n  =  0: 


a®„  n  _i_  a£»  ,  £^of 

dn  “  r  (3)  da  ^  dr] 


■■  r  (3)  os. 


a®,  *** 


'0, 


(1.25) 


a3  dr\ 

Prom  the  first  relationship  of  Conditions  (1.25)  we  see  that  the 
standard  derivative  of  the  harmonic  function  $>q  is  zero  at  the  bound¬ 
ary  of  the  cavity.  Hence  it  is  possible  to  consider  O  =  const  without 
.  limiting  the  generality. 

The  components  of  the  vector  1  in  the  curvilinear  system  of 
coordinates  are  associated  with  the  components  of  the  same  vector  in  , 
the  Cartesian  coordinate  system  by  Relationships  (1.19).  These  rela¬ 
tionships  are  linear  with  regard  to  ¥q2,  f and  the  coefficients 

of  Yq2,  Yqj  are  Independent  of  T]j  consequently  the  functions  ¥qq, 

and  TQn  satisfy  the  same  equations  a3  YQ1.  Thus  : 

a*rv 

®'^0|S  =  > 


gtry 

-w  _  To« 
- ~frf  ' 


^  ■  oT—S  •  (1.2^) 


dr 


At  the  boundary  of  the  cavity  these  functions  satisfy  the  follow¬ 
ing  boundary  conditions: 


dr] 


■  =»  0, 


(1.25) 


The  general  solution  of  the  equation  satisfied  by  Tq^  and 


YQn  has  the  following  form: 


«  -  «/*<+  e*-*' 


Let  Re  /a"*  >0.  By  assumption,  far  from  the  walls  of  the  cavity 
there  are  no  vortices.  Just  as  is  done  in  boundary -layer  theory,  we 
will  consider  the  value  tj  =  «  to  be  the  respective  internal  points 
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of  the  cavity  sufficiently  far  from  the  boundary.  Then,  In  agreement 
with  the  assumption  concerning  the  absence  of  vortices,  for  -  «#  we 
will  have  c^  =  0. 

The  arbitrary  constant  of  Integration  Cg  is  determined  from 
/Conditions  (1.25).  Finally,  for  the  functions  ^On  we  obtaln: 

•'Vo.-'O.  '  -  -  r  (P)  a  =  0  (1.26) 

We  will  limit  ourselves  to  determination  of  only  the  first  terms  of 
Series  (1.21).  Within  the  limits  of  this  accuracy  we  write  the 
components  of  the  absolute  velocity  vector  V'  in  a  stationary  system 
of  coordinates  with  its  origin  at  the  point  of  suspension  of  the 
pendulums 

•  u'  =  c-5-  [/a  —  r  (p)  a  cosa  e~Va1i  -f  O  (e)]  e"‘ 

1  _  '  w'  =  0  .  . 

v'  =-c|-  [r  (P)  o  sina  e~y°ri  +  0  (e)]  e8^  .(1*27 ) 

where  l  is  the  distance  of  the  center  of  gravity  of  the  cavity  from 
the  suspension  axis,  relative  to  the  characteristic  dimension  of  the 
cavity.  The  resulting  solution  (1.27)  will  be  asymptotic.  Relative 
to  it,  we  have  the  following  result,  proof  of  which  will  not  be  cited 
here.  The  difference  modulus 

lv'-Vo'|<0(e) 

Here  Y*  is  the  exact  solution  of  the  linearized  problem  and  Y*0 
is  the  approximate  solution  of  this  problem  obtained  by  us. 

Thus,  if  the  Reynolds  number  is  sufficiently  high,  the  approxi¬ 
mate  solution  assures  good  accuracy. 

2.  For  a  final  solution  of  the  problem  it  is  necessary  to 
determine  the  as  yet  unknown  constant  a.  For  this  it  is  necessary 
to  use  the  equation  of  pendulum  oscillations  which  may  be  obtained 
using  the  theorem  concerning  the  variation  of  angular  momentum  of  a 
system.  We  have: 
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(2.1) 


j};\  (r'  x  PiVj')  *  +  Jj  (  (r'  x  pV')  dx  =  (  (r'  x  ?&)  dx  +  ^  (r'  x  pg)  dx 
d,  8  8,  d 


Here  -  volume  of  the  solid; 

D  =  volume  of  cavity; 

and  T[  =  density  and  velocity  of  points  of  the  body; 

p  and  y  =  density  and  velocity  of  a  particle  of  fluid; 

g  =  vector  of  gravitational  acceleration; 

r!  =  radius  vector  from  the  axis  of  rotation  to  the 
point  of  the  body  or  fluid. 

In  our  case  the  equation  of  moments  (2.1)  yields  a  projection, 
which  differs  from  the  identical  zero,  only  onto  the  z1  axis. 

The  integrals  which  enter  into  (2.1)  [returning  to  dimensional 
variables  according  to  Formula  (1.5)]  may  be  reduced  to  the  following 
form 

V"  (p-f)  (2-2) 

|fL  J  (r-  x  pV')  rft}x.  -  c  (Mi' V*  +  V^  MQ o' 1 )  e"*1'  [mQ  -  2nP  J  r«  (P  W)  (  2 . 5  ) 

{ ^  (*'  *  Pig)  dx  +  ^  (r'  x  pg)  dt|z,  =  —  eg  ( Mtlx  +  Ml’)  e*'r  (2.4) 


Here  *=  mass  of  the  solid; 

M  =  mass  of  the  fluid; 

k  =  radius  of  inertia  of  the  solid  relative  to  the  axis  of 
suspension  of  the  pendulum; 

l *  =  distance  from  the  axis  of  suspension  to  the  center  of 
mass  of  the  cavity; 

IJ  =  distance  from  the  axis  of  suspension  to  the  center  of 
gravity  of  the  solid; 

(3£,  «  coordinates  of  the  poles  of  the  cavity. 

Substituting  (2.2),  (2.5),  and  (2.4)  into  (2.1)  we  obtain  the 
following  equation  to  define  a; 

3'«  +  1 1/7 o' 1 '  +  a>«  =■  0  •  (x  -  ’“'-r  S%\  +  M*)  (2-5) 

A 

We  will  set  a  *  vX  x.  Then  (2.5)  takes  the  form: 
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X*  +  3*  +  JijL  =  0  ^han  y*  +  y3  +  q*  ==  0  (Vx  -  tjr,  »  ?J)  (2.6) 

We  will  prove  that  this  equation  has  only  two  roots  which  satisfy 
the  condition  Re  y  >  0.  We  set  y  »  a  +  ifi; 

After  substitution  in  (2.6),  we  obtain  for  the  determination  of 
a  and  0  the  system 

(«*  “  PV  -  4a’P*  +  a*  -  3a?»  +  q'  =  0  ,  -  _  v 

, .  4a>  -  4aP‘  -  p*  +  3a*  «  0  '  ' 

p 

Determining  £  from  the  second  equation  of  System  (2.7)  and 
substituting  it  in  the  first  equation  we  obtain: 

64a*  +  96a‘ +  48a*  +  8a3  -  9*  (16d‘ +  8a  +  1)  =  0.  ^  %  -  ±  a j, /~|» ± |  ('2.8) 

The  first  equation  in  (2.8)  has  only  one  change  of  sign  of  the 
coefficients.  According  to  Descart's  rule  this  polynomial  has  one 
positive  root.  We  are  not  interested  in  negative  roots  since  in  this 
case  the  condition  Re  /x~  >  0  is  violated.  The  second  equation  of 
(2.8)  gives  two  values  of  p  which  differ  only  in  sign.  Therefore, 

Eq.  (2.6)  has  only  two  roots  for  the  case  Re  y  >  0. 

We  will  solve  Eq.  (2.6)  for  the  two  limiting  cases:  when  £ 
is  large  and  when  £  is  small.  Let  us  recall  that  the  solution  of 
the  problem  concerning  the  motion  of  a  fluid  within  a  pendulum  cavity 
was  obtained  for  a  large  Reynolds  number.  In  turn,  the  magnitude  of 
the  Reynolds  number  depends  essentially  on  the  characteristic  period 
of  one  oscillation  tr/  j  Im  o7  |  which  is  determined  only  after  solu¬ 
tion  of  Eq.  (2.6).  Consequently,  of  all  the  possible  solutions  of  Eq. 
(2.6)  which  correspond  to  different  parameters  of  the  pendulum  it  is 
possible  to  use  only  those  which  guarantee  a  sufficiently  large 
Reynolds  number.  Later  we  will  see  that  in  the  limiting  cases  of 
large  and  small  £  which  we  examined,  it  is  always  possible  to  select 
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pendulum  parameters  such  that  the  Reynolds  number  remains  large. 

Let  us  examine  the  case  of  large  £.  We  will  seek  the  solution 
of  Eq.  (2.6)  In  the  form  of  a  series  In  powers  of  1/ /"q" 

V-Vl{lh  +  y*V  !  +  •••)  (2.9) 

Substituting  (2.9)  in  (2.6)  and  equating  to  zero  the  sum  of  the 
coefficients  with  identical  powers  of  l/  /q~  we  obtain  equations, 
the  determination  of  the  unknowns  y0,  y^,  etc. 

For  y0  we  haves 

y»*  +  1  -  0,  Rey0>0 

whence  we  obtain: 

=  ~z~  ^  0*  **  "j-  (l  —  0 

For  the  determination  of  y^  we  have 

4y,  +  1  «*  0  when  yx  —  —  V« 

With  an  accuracy  to  value  of  the  order  of  l/  /q-  inclusive,  we 


find 


or 


y 


VZ 

~T 


FI#) 

(o'  =»  vi*y*) 


(2.10) 


Thus,  in  the  case  under  consideration  the  frequency  n  and  ampli¬ 
tude  A  of  the  pendulum  oscillations  will  be,  respectively. 

As  v->  0,  the  damping  factor  approaches  zero  while  the  oscillation 
frequency  approaches  the  frequency  co  which  corresponds  to  the  oscllla-  - 
tion  of  a  pendulum  with  an  ideal  fluid.  We  will  calculate  by  what 
amount  the  amplitude  decreases  during  each  period  T  ■  ir/  |lm  o'  |. 
Denoting  this  amount  by  A  we  have. 
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A  =  1 


exp  Re  (a 'h')  _  n  Res' 

exp  Re  (o'*!')  “  v  llmo'l 


(2.11) 


Determining  Re  a'  and  |lm  a*  |  from  (2.10)  and  substituting  In 


(2.11)  we  obtain: 


A  -  1  — exp  (-aTJpr)*^  (A -.0  when  ?-<») 

It  Is  easy  to  see  that  the  examined  case  of  large  £  corresponds 
to  oscillations  close  to  those  of  a  pendulum  with  an  ideal  fluid. 

Let  us  consider  the  simplest  example.  Let  the  pendulum  be  a 
weightless  spherical  shell  filled  with  a  viscous  fluid  and  connected 
to  the  suspension  axis  by  a  weightless  rod  (Fig.  3).  For  such  a  pendulum 
we  will  have: 


(2.12) 


The  characteristic  period  of  one  oscillation  with  an  accuracy 
to  values  of  the  order  of  l/  /q  Is 


T-^VTfJ  (2.13) 

The  magnitude  of  T  Increases  with  an  Increase  In  the  length  l 

of  the  pendulum.  In  this  case  there  Is  the  danger  of  an  Inadmissible 

reduction  of  the  Reynolds  number.  We  will  set  the  Reynolds  number 

4 

%e  an*  evaluate  the  allowable  pendulum  lengths  l  for  various 

values  of  radius  R.  By  definition  NRe  =  R2/vT. 

4 

In  order  to  keep  NRe  ^  10  It  Is  required  "that 


T**  RX  ^  1  R*t 

^vlO*  0r 

Suppose  that  the  pendulum  is  filled  with  water  at  a  temperature 
of  20°C  (v  «■  1.01*10“^  m2/sec).  For  R  =  0.1  m  we  have  l/R  £  10, 
i.e.,  I  £  1  m.  For  R  -  1  m  we  obtain  l/R  <  lo\  i.e.,  i  £  10  km. 
With  an  Increase  in  R  the  allowable  value  of  l/R  increases  as  R? . 
Thus  the  range  of  applicability  of  solution  (2.12)  is  sufficiently 

“  .  *  . I  ‘  ■  *■.*'-*■ . V  “ » 
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broad. 


Let  us  now  examine  the  case  of  extremely  small  £.  We  note  thr 
it  is  possible  to  realize  small  £  not  only  by  increasing  the  kinen: 
viscosity  v,  which  might  lead  to  a  reduction  in  NRq.  Actually,  in 
the  case  of  the  pendulum  Just  considered  it  is  apparent  that  this  n.«. 
be  achieved  by  decreasing  the  length  of  the  pendulum  l  or  Increasing 

< 

the  radius  of  the  cavity  R. 

Let  us  represent  Eq.  (2.6)  in  the  following  form: 


y'  (y  + 1)  -  -  9* 


(2.14) 


From  this  it  follows  that  , 

|y*l|y  +  l|  =  l- ?*l=  ?* 

The  product  |y^|  |y+l j  decreases  with  decreasing  In  this 

process  the  factor  |y+l|  cannot  approach  zero,  otherwise  the  condition 

Re  y  >  0  is  violated.  Consequently  the  modulus  of  the  unknown  root 

2/"3 

has  an  order  of  smallness  q  '  y .  Starting  from  this  we  will  seek  the 
solution  of  Eq.  (2.6)  in  the  form  of  the  following  series: 


y  —  9V*  (l'o  +  q'tyi  +  . . .) 


For  the  determination  of  the  unknowns  yQ  and  y^  we  obtain 

y0*  +  1  =  0,  y*  4-  dyl  =  0 


From  thiB  we  find: 

J/01  =  v  +  ‘  v  ■  yw 


v — t  x  ’  y« 


-  A  1  A  » 


yn 5 


i+ilS' 

5  '  « 


,2/5 


With  an  accuracy  to  values  of  the  order  q  inclusive,  we  have: 

v  =  7  ?7‘  K‘  +  7  «*)  ±i/3(i-i  ?v.)j 


For  a'  we  obtain  the  expression: 

o>  =  -E.- r_f  1_  ±1^31 

2vV*l  \  3  }±l  —  \ 

or,  limiting  ourselves  to  magnitudes  of  the  order  q 


(2.15) 


V?. 


Just  as  in  the  preceding  case,  we  determine  by  what  amount  the 
amplitude  drops  from  its  original,  value  during  the  characteristic 
period  T. 


FTD-TT-65-5H/I+2+4 
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We  have 


A  =»  1  —  exp  n  Re  o'  /  |  Im  «'  | 

Using  Eq.  (2.15)  to  determine  Re  a1  and  llm  o*  j,  we  obtain 


^A— *•! — exp( — ^=j»0.84  wh«n  q  —+0j 


(2.16) 


It  is  apparent  that  this  case  differs  greatly  from  the  case  of 
oscillations  of  a  pendulum  ‘with  an  ideal  fluid.  It  turns  out  that  it 
is  possible  to  select  a  pendulum  cavity  so  large  or  a  pendulum  length 
so  small  that  the  influence  of  viscosity  on  the  motion  of  the  pendulum 
cannot  be  neglected  despite  large  Reynolds  numbers.  In  the  case  of 
extremely  small  £  the  amplitude  may  drop  on  the  order  of  84#  from  its 
original  value  during  a  single  swing. 

We  will  demonstrate,  using  as  our  example  the  pendulum  shown  in 
Pig.  J>,  that  by  increasing  the  radius  of  the  cavity  R  while  keeping 
l  and  v  constant  it  is  possible  to  reach  as  large  NRe  and  as  small 
£  as  desired.  For  this  pendulum  we  have: 

(*•»> 

4W  (2.18) 


g-VjILJL, 

*  4v  ST' 


N  -^*V\ 

Nr*  "ST ' 


It  is  easily  seen  from  (2.18)  that  NRe-»  w  and  q-»  0  as  R -*  «. 
Oscillations  of  the  pendulum  in  this  case  are  extremely  different 
from  oscillations  of  a  similar  pendulum  with  an  ideal  fluid.  Thus, 
a  pendulum  filled  with  water  at  20°C  (v  -  1.01 *10“^  m2/sec)  with  a 
cavity  radius  R  »  0.1  m  and  a  length  l  -  0.00016  m  will  complete  5 
oscillations  per  second,  while  the  same  pendulum  with  an  ideal  fluid 
will  make  25  oscillations  per  second.  For  the  pendulum  parameters 
indicated  above  we  will  have: 
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FTD-TT-63-511/1+2+4 


?  as  0.0012,  NRt  as  10* 

which  completely  justifies  the  use  of  Solution  (2.17)  in  this  case. 


Submitted  January  29,  1962 


Pig.  1.  Fig.  2.  Fig.  5. 
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